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ABSTRACT 

Multi- frequency  modulation  is  a  highly  bandwidth  efficient 
signalling  technique  for  digital  communications.  In  order  to  make 
the  technique  as  insensitive  as  possible  to  unknown  or  fluctuating 
phase  and  amplitude  changes  in  the  channel  transfer  function 
between  the  transmitter  and  receiver,  frequency  domain  differential 
encoding  techniques  have  been  developed  and  their  error  performance 
calculated.  It  is  shown  that  in  the  range  of  2  to  8  bits/Hz  of 
channel  bandwidth  efficiencies,  differential  encoding  results  in  a 
penalty  of  3  to  5  db  in  required  E^Nq  when  compared  to  fully 
coherent  multi-frequency  modulation.  Design  procedures  are 
presented  that  provide  near  optimum  QAM  constellations  for  fully 
differential  coding  and  for  the  hybrid  scheme  of  differential  phase 
coding  and  absolute  coding  of  amplitudes. 
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I  INTRODUCTION 

Multi-frequency  Modulation 

Multi-frequency  modulation  (MFM)  is  a  digital  modulation 
technique  that  is  being  researched  at  the  Naval  Postgraduate  School 
for  application  in  military  digital  communications  systems.  It  is 
a  high-data-rate,  highly-bandwidth-efficient,  modulation  method 
with  important  potential  for  point-to-point,  broadcast  and  packet 
switched  network  digital  communications  links.  The  method  makes 
extensive  use  of  digital  signal  processing  techniques  for  its 
implementation . 

MFM  is  also  known  as  orthogonal  frequency  division 
multiplexing  (OFDM) ,  because  of  the  orthogonality  of  the 
multiplicity  of  carrier  frequencies  that  are  used  to  transport  the 
digital  data.  OFDM  has  been  discussed  in  the  literature  in  the 
context  of  data  modems  for  the  public  switched  telephone  network 
(PSTN)  ,  as  a  groupband  modem,  as  a  mobile  HF  modem,  and  for 
satellite  broadcast  of  digital  audio  programs (DAB) [1],[2],[3],[4]. 
One  of  the  difficulties  in  the  practical  application  of  MFM  (OFDM) , 
is  that  of  carrier  recovery  for  demodulation.  Unlike  ordinary 
bandpass  digital  modulation  techniques  that  employ  a  single 
carrier,  MFM  uses  a  great  number,  perhaps  hundreds  or  thousands,  of 
carriers  simultaneously.  If  the  channel  contains  any  envelope 
delay,  that  is  phase  shift  that  is  not  linear  with  frequency,  then 
the  channel  transfer  function  must  be  known,  or  estimated, 
accurately  if  the  carrier  phases  are  to  be  used  to  carry  data. 


During  the  past  year  we  have  concentrated  on  solving  this 
problem.  A  solution  that  we  have  found  works  well  in  a  practical 
channel  [5],  is  to  differentially  encode  information  in  the  phase 
of  adjacent  tones  in  the  same  baud.  The  purpose  of  this  report  is 
to  document  the  theoretical  performance  of  this  and  of  more  general 
classes  of  MFM  differential  modulation. 

The  underlying  principles  of  MFM  generation  and  demodulation 
were  reported  last  year  [6].  The  results  of  that  report  serve  as  a 
point  of  departure  for  this  work.  An  MFM  packet  is  shown  in  Figure 
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Figure  1  An  MFM  Packet 
1.  The  terms  used  in  this  report  are  consistent  with  those  of  [6] 


Siqnal-to-Noise  Ratio  and  Eb/N0 

It  was  shown  in  our  previous  report  [6,  pg  24]  that  for  any 
given  baud  corrupted  by  additive  white  gaussian  noise (AWGN) ,  the 
output  of  the  DFT  in  the  MFM  receiver  is  the  complex  frequency 
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domain  sequence 

Y(k)=  S(k)  +  W(k)        ;  Jc,  <  k  <  k2         (l.l) 

where  the  received  signal  baud  is  given  by 

S(k)  =   Akcos0k   +jAksin0k  ,  (1.2) 

and 

Ak  =kx(2Pk)V2 
with  Pk  the  received  tone  power  and  kx  is  the  number  of  points  in 
the  DFT.   W(k)  is  a  white  zero  mean  complex  gaussian  frequency 
domain  sequence  with  uncorrelated  real  and  imaginary  parts  that 
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have  equal  variances  ^AfNj/4.  The  frequency  spacing  between  tones 
is  Af.  By  dividing  (1.1)  by  the  standard  deviation  of  the  noise 
components  one  obtains  the  unit  variance  complex  sequence 

F(k)  =  R(k)  +jl(k)  =  S'(k)  +  W'(k)  (1.3) 

with  mean  value 

E[F(k)]    =   S'(k)    =        A'kcos0k  +jA'ksin0k  (1.4) 

where 

A\   =(2Ek/N0)*    .  (1.5) 

We  have  made  use  of  the  fact  that  Pk/Af=Ek  is  the  received  energy 
in  the  kth  tone.  This  is  because  Af  is  the  reciprocal  of  the  baud 
length  AT. 


In  some  applications,  the  channel  signal-to-noise  ratio  is  the 
parameter  most  pertinent  to  system  performance.  We  shall  refer  to 
the  signal-to-noise  ratio  in  the  signal  band  as  the  wideband 
signal-to-noise  ratio,  SNRyB,  and  define  it  as  the  average  signal 
power  in  the  band  to  the  average  noise  power  in  the  same  band.  Thus 

SNRyB  =  KP/N0W  =  E/N0  (1.6) 

where  P  and  E  are  the  average  tone  power  and  tone  energy 
respectively,  and  K  is  the  number  of  tones  in  the  signal  band  of 
width  W  =KAf .  Furthermore,  if  there  are  Q  bits  transmitted  per  tone 
then   Eb  =  E/Q  is  the  average  energy  per  bit  transmitted  and 

Eb/N0  =  SNR^/Q  =  E/QN0.  (1.7) 

We  shall  use  (1.5)  through  (1.7)  freely  when  making  performance 
comparisons  of  various  types  of  MFM  modulation. 

MFM  Constellations 

Three  classes  of  modulation  will  be  considered  in  this  report; 
phase  modulation  (PSK) ,  amplitude  modulation  (ASK)  and  combined 
amplitude  and  phase  modulation.  We  shall  refer  to  the  latter, 
hopefully  unambiguously,  as  QAM.  (Ambiguity  could  arise  from  the 
fact  the  the  PSK  is  coded  and  decoded  using  its  guadrature 
components,  but  the  context  of  our  discussions  should  make  the 
distinction  between  QAM  and  straight  PSK  clear.)  In  each  case  we 


shall  present  results  for  both  differential  and  non-differential 
encoding.  A  particuarily  useful  digital  modulation  format  is  a 
hybrid:  differential  phase  coding  and  non-differential  amplitude 
coding.  The  hybrid  constellation  for  MFM  has  been  developed  by 
Mercury  Digital  Communications,  Inc. [7]. 

Differential  modulation  is  obtained  by  encoding  the  input  data 
symbols  as  changes  in  signal  state  occupancy  of  adjacent  tone 
constellations  in  such  a  way  that  the  changes  are  invariant  under 
a  rotation  of  the  constellations.  The  principal  purpose  for 
differential  encoding  is  to  render  the  transmission  insensitive  to 
uncompensated  channel  phase  shift.  In  MFM  the  tones  can  be  placed 
very  close  together.  If  adjacent  tones  are  close  enough  together 
then  their  uncompensated  phase,  that  is  the  rotation  of  their 
constellations,  will  be  nearly  identical,  and  it  cancels  out  in  the 
decoding. 


II  ERROR  PATTERNS  AND  PERFORMANCE  MEASURES 

Symbols  and  Bits 

In  MFM,  Q  bit  symbols  are  coded  into  the  phase  and  amplitude 
of  each  of  the  multiplicity  of  carrier  frequency  tones  during  each 
baud.  The  symbols  are  partitioned  into  M  bits  that  are  coded  into 
tone  phase  and  N  bits  that  are  coded  into  tone  amplitude  (Q=M+N) . 
We  adopt  the  notation 

s  =  p_,a  (2. 1) 

to  designate  the  symbol  and  its  partition  into  phase  and  amplitude 
bits.  For  example,  with  4  bits  per  tone,  one  may  have 

1010  =  101,0  (2.2) 

and  three  bits  of  phase,  in  this  example  101,  are  coded  as  one  of 
eight  possible  tone  phases  and  one  bit  of  amplitude,  in  this 
example  0,   is  coded  as  one  of  two  possible  tone  amplitudes. 

Errors  may  occur  in  the  receiver  in  any  of  the  Q  bit 
positions.  An  error  pattern  is  defined  as  a  Q  bit  symbol 


e  =  ep,ea  (2.3) 


which  we  also  partition  into  phase  bit  errors  and  amplitude  bit 

Q  M 

errors.  There  are  2  total  error  patterns  made  up  from  2  phase 


u 

error  patterns  for  each   of  the  2  amplitude  error  patterns 
The  output  of  the  receiver  is  the  decoded  symbol 


d  =  s  +  e  (2.4) 

where  the  plus  sign  designates  bit  by  bit  binary  addition  (logical 
"or")  of  the  transmit  symbol  with  the  error  pattern.  The  decoded 
symbol  is  also  partitioned  into  decoded  phase  and  amplitude  bits, 


d  =  dp,da  (2.5) 


A  symbol  error  occurs  anytime  that  d  *  s,  that  is  for  any  of  the 
2  -1  non-zero  error  patterns,  e  *    0. 


Symbol  Error  Probabilities  and  Bit  Error  Rate  (BER) 
The  probability  of  symbol  error  is 

Ps  =  2  Pr[e#0|s=si]Pr[s=si]   .  (2.6) 


.  Q 

If  we  assume  proper  source  coding  then  each  of  the  2  data  symbols 


will  be  equally  probable  and 


Ps  =  2"°  S  Pr[e^o|s=si]  .  (2.7) 

i 


Performance  of  digital  communications  systems  is  customarily 


specified  is  terms  of  bit  error  rate  (BER) .  BER  is  defined  as  the 
expected  number  of  bit  errors  in  a  very  long  sequence  of  symbols 
divided  by  the  number  of  bits  in  the  sequence.  For  example,  if  the 
expected  number  of  bit  errors  in  a  sequence  of  one  million  bits  is 
one,  then  the  BER  is  10  .  As  long  as  symbols  that  make  up  a 
sequence  of  bits  are  statistically  independent  and  the  error 
probabilities  are  stationary,  then  the  BER  is  also  the  expected 
number  of  bit  errors  in  any  given  received  symbol  divided  by  Q,  the 
number  of  bits  in  the  symbol.  We  shall  adopt  this  definition  of  BER 
in  this  report. 

Let  B  be  the  number  of  symbol  bits  in  error,  that  is  the 
number  of  ones  in  e.  Then 

E[B]  =  2  E[B|s=s1-]Pr[s=si]  (2.8) 

and  again  assuming  equally  probable  symbols, 

E[B]  =  2"°  £  ECBls^sJ  (2.9) 

with 

BER  =  E[B]/Q  (2.10) 

in  either  case. 

Single  and  Multiple  Bit  Errors 

BER  is  a  normalized  expected  value.  Sometimes  it  is  most 
convienent  to  compute  it  as  a  weighted  superposition  of  terms 
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corresponding  to  single  bit  errors,  double  bit  errors,  etc,  up  to 
Q  bit  errors.  That  is, 


Q  -a    *         . 

E[B]  =  2  b  Pr[B=b]  =  2   S  (2  b  Pr  [  B=b  |  s=s,.  ]  }    (2.11) 


where,  Pr[B=b|s=si]  are  the  conditional  bit  error  probabilities. 
These   can   be   computed   from   the   conditional   error   pattern 
probabilities  as  follows; 

Pr[B=b|s=si]  =   2  Pr[e|e  contains  b  oneSjS^-]    (2.12) 

where  the  sum  is  over  the  Nb  =  (Q,b)  error  patterns  containing  b 
ones  out  of  Q  bits.  For  example,  there  are  N1  =  (Q,l)  =  Q  single 
bit  error  patterns,  N2  =  (Q,2)  =  Q(Q-l)/2  double  bit  error  patterns 
etc. 

In  summary,  the  symbol  error  probability  is  computed  from 
(2.7)  and  the  BER  from  (2 . 10) , (2 . 11)  and  (2.12).  Both  computations 
require  the  conditional  error  pattern  probabilities  Pr[e=e}  |  s=si  ]  . 
The  objective  of  the  next  three  Chapters  is  to  obtain  theoretical 
expressions  for  these  for  various  MFM  and  MFDM  constellations. 


Ill  PHASE  MODULATION 

In  this  Chapter  MFM  and  MFDM  systems  that  use  constant 
amplitudes  and  employ  only  phase  modulation  to  encode  the  Q  symbol 
bits  (M=Q,  N=0)  will  be  considered.  Two  types  of  modulation  will  be 
analyzed;  coherent  modulation/demodulation  and  differential 
modulation/demodulation.  In  coherent  modulation,  the  absolute  phase 
of  each  of  the  carrier  tones  must  be  known  at  the  receiver.  This 
form  of  modulation  should  only  be  used  in  a  linear  phase  channel, 
or  in  a  channel  that  has  been  equalized  to  linear  phase.  Any 
residual  phase  jitter  will  affect  performance  as  noise,  and  any 
significant  unequalized  non-linear  phase  (group  delay)  will  most 
probably  render  coherent  phase  modualtion  unworkable  at  any  SNR. 
However,  its  analysis,  assuming  perfect  phase  equalization,  is 
included  for  comparison  with  differential  modulation. 

Differential  phase  modultion  of  MFM,  MFDM,  encodes  phase  bits 
as  the  phase  difference  between  adjacent  carrier  frequency  tones  in 
the  same  baud.  Although  it  is  not  necessary  that  adjacent  tones  be 
used,  or  that  the  same  baud  be  used  (we  could  encode  using  the 
phase  difference  of  each  of  the  same  tones  on  successive  bauds  as 
is  done  in  conventional  single  carrier  differential  phase 
modulation  (DPSK)  transmission) ,  the  theory  behind  differential 
encoding  in  the  frequency  domain  using  adjacent  tones  is  that  the 
differential  phase  jitter  between  adjacent  tones  introduced  by  the 
channel  can  be  reduced  to  any  desired  level  by  putting  the  tones 
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sufficiently  close  together  in  frequency.  Recall  that  the  tone 
spacing  in  MFM  is  Af  =  1/AT.  Therefore  differential  phase  shift 
through  an  arbitrary  unequalized  channel  is  inversely  proportional 
to  baud  length. 

Coherent  Phase  Modulation/Demodulation 

M 

A  phase  constellation  for  2  -ary  PSK  is  generated  by  utilizing 
a  constant  tone  amplitude  and  2  phases.  The  phases  are  equally 
spaced  at  intervals  A0=  27T/2  ,  then  rotated  by  adding  A0/2.  Thus 
the  modulation  phase  bits  are  encoded  according  to  the 
transformation, 


pfk)-  -  d0(k)  =  i(k)A0    ;   i  =   0,1, 2M-1 


(3.1) 
0(k)  =  d0(k)  +  A0/2 


Such  a  constellation  for  8  phases  is  shown  in  Figure  2. 

The  maximum  likelihood  demodulation  rule  (MLD)  is  to  decode 
the  complex  point  Y(k)  generated  by  the  receiver  DFT(the  matched 
filter  receiver  for  MFM,  see  [6])  to  the  nearest  point  of  the 
constellation.  The  decoding  boundaries  for  MLD  are  sectors  of 
anglular  width  A0  as  illustrated  for  an  2  of  eight  in  Figure  2.  We 
shall  refer  to  this  type  of  decoding  as  "sector  decoding"  (SD)  . 
Notice  that  SD  is  independent  of  tone  amplitude  A,  or  equivalently 
of  E/N0.  This  makes  the  implementation  of  SD  in  the  receiver 
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independent  of  system  gain,  a  highly  desireable  feature.  Of  course 
system  performance  is  not  independent  of  SNR. 


i 
X  0\;1 

X^^       • 

00  1      / 

,/        •  oo,o 

11,0    •           / 

Figure  2   An   Eight   PSK   Constellation  with   Sector   Decoding 
Boundaries 


a.  Symbol  Error  Probability 

Phase  symbol  error  decoding  probabilities  are  the  same  for 

M 

each  of  the  2  symbols.  This  is  so  because  the  statistical 
properties  of  the  complex  noise  W(k)  are  invariant  to  phase 
rotations.  That  is,  the  real  and  imaginary  parts  of  W(k)exp(j0) 
are,  like  those  of  W(k) ,  uncorrelated,  equal  variance  gaussian 
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random  variables.  Consequently,  the  probability  of  symbol  decoding 
error  given  by  (2.7)  is  the  same  as  the  conditional  probability  of 
decoding  error 

Ps  =  Pr[e#o|  s^.]  ,  for  any  i.  (3.2) 

The  probability  that  e  *  0  when  s  =  s,-  is  the  probability  that 
Y(k),  or  equivalently  F(k) ,  is  not  in  the  sector  of  angular  width 
A0  defined  by  the  angles  iA0  and  (i+l)A0. 

i.)   QPSK  (M  =  2) 

We  include  this  case  for  completeness  as  the  results  are 
well  known  and  were  presented  in  the  previous  report [6].  Also,  for 
more  phase  states  and  narrower  phase  sectors,  the  results  that  will 
be  presented  in  the  next  section  are  only  approximate (the 
approximation  becomes  more  accurate  as  the  sectors  become  more 
narrow) .  Exact  results  are  easily  found  for  QPSK  as  follows: 

The  probability  that  F(k)  is  not  in  the  first  quadrant, 
given  s^  is  simply  one  minus  the  probability  that  F(k)  is  in  the 
first  quadrant.  The  probability  that  F(k)  is  in  the  first  quadrant 
is  the  probability  that  R(k)>0  and  that  l(k)>0.  Since  the  real  and 
imaginary  parts  of  F(k)  are  statistically  independent; 

Ps  =  1  -  Pr{R(k)>0}Pr{I(k)>0} .       (3.3) 
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These  probabilities  are  easily  found  using  the  fact  that  R(k)  and 
I(k)  are  unit  variance  gaussian  random  variables  with  mean  values 
determined  from  (1.4) -(1.6)  to  be 


with 


E{R(k)}  =  E{I(k)}  =  A\/(2)*  "  P  (3-4) 


p  =  {SNI^}*  =<E/N0}* 


where  we  have  used  the  fact  that  <pk   =   it /A    for  s  =  Sq,  the  symbol 
assumed  as  transmitted   on  tone  k.  Consequently, 


Ps  =  1  -[1-  Q(p)]2   «  2Q(p)  (3.5) 


is  the  symbol  error  probability.  The  approximation  is  very  good  for 
SNRyB  greater  than  about  10  db.  Note  from  (1.7)  that  for  QPSK 


p2  =  E/N0  =  SNRyg  =  2Et/No-  (3-6) 


ii. )  2H-ary  PSK 

Let  p(A9)  be  the  probability  that  F(k)  falls  outside 
an  angular  sector  of  width  ±Ae/2  with  respect  to  0i ,  the  angle  of 
F(k)  given  the  transmission  of  s,-  in  the  absence  of  noise.  That  is, 

p(A6)  =  Pr[|o  -  0,|>A6/2]  (3.7) 


14 


where  a   is  the  angle  of  Y(k)  (or  equivalently  of  F(k)). 
It  is  known  that  for  Ae  <  n/2,      this  probability  is  well 
approximated  by  ([8];  Hakin  pg  317). 


p(A9)  «  2Q{A'ksin(Ae/2)  }  =   2Q{  (2)*psin(Ae/2)  }   (3.8) 


and  that  the  approximation  becomes  increasingly  more  accurate  with 
increasing  p  and  decreasing  AG.  A  PSK  symbol  error  occurs  if 
|  a  -  0i  | >A0/2 .  Therefore,  the  PSK  symbol  error  probability  is 


Ps  =  p(A0)  m  2Q{(2)Vin(A0/2)}.  (3.9) 


We  note  that  (3.9)  with  A0  =  7r/2,  which  is  QPSK,  agrees  with  the 
approximation  generated  by  binomial  expansion  of  the  exact  equation 
derived  for  QPSK  and  presented  in  (3.5). 


Differential  Phase  Modulation/Demodulation 

In  MFDPSK,  the  phase  bits  for  tone  k  are  encoded  into  equally 
spaced  phases 


p_ikli  -*  d0(k)  =  i(k)A0    ;  1=1,2, 2M-1,     (3.10) 


and  this  differential  phase  shift  is  added  to  the  phase  of  the 
previous  tone  to  obtain  the  tone  phase  of  tone  k  as 
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0(k)  =  0(k-l)  +  d0(k).  (3.11) 

An  initial  condition  must  be  chosen.  By  choosing  the  phase  of  the 
first  carrier  frequency  tone  as 

0(k1)  =  A0/2  (3.12) 

the  constellations  of  all  tones  are  rotated  by  A0/2  thereby 
yeilding  transmit  tone  phases  identical  to  those  used  in  coherent 
MFPSK  discussed  in  the  previous  section. 

Differential  decoding  the  phase  bits  is  accomplished  by 
computing  the  complex  rotated  product 

D(k)  =  F(k)F*(k-l)exp[jA</>/2]  ;      k^l^fc^kj  (3.13) 

for  adjacent  digital  frequencies  in  each  MFDM  baud.  In  the  absence 
of  noise,  the  phase  of  D(k)  is  the  phase  difference  between  tone  k 
and  tone  k-1  plus  A0/2,  that  is  i(k)A0  +  A0/2.  Comparison  with 
(3.1)  reveals  that  the  complex  rotated  products  have  exactly  the 
same  constellations  as  described  for  coherent  PSK.  The  differential 
receiver  optimally  decodes  the  phase  bits  at  this  stage  by  using 
the  same  sector  decoding  rules  described  above, 
a.  Statistical  properties  of  the  Complex  Rotated  Product. 
Define  the  complex  random  product 


Fp(k)  =  F(k)F  (k-1)  =  S'(k)S'  (k-1)  +  Wp(k)    (3.14) 
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with 


Wp(k)  =  W  (k)W'*(k-l)  +  S»  (k)W*(k-l)  +  W  (k)S»*(k-l)  .   (3.15) 


Since  W  (k)  is  zero  mean,  it  follows  directly  that  W  (k)  has  zero 
mean  and  therefore  that  the  mean  of  Fp(k)  is 


E[F_(k)]  =  S'(k)S'*(k-l)  =  {2(EkEk.1)VNo}exp[jd0(k)].      (3.16) 


The  real  and  imaginary  parts  of  the  complex  product  are 


Rp(k)=  Re[F(k)F  (k-1) ]  =  R(k)R(k-l)  +  I(k)I(k-l) 


Ip(k)=  Im[F(k)F  (k-1)]  =  I(k)R(k-l)  -  R(k)I(k-l). 


(3.17) 


Because  of  the  statistical  independence  between  adjacent  tones  and 
between  the  real  and  imaginary  parts  of  the  same  tone,  it  follows 
directly  that 

Var{Re[Wp(k)  ]  }=  Var{Rp(k)  }  =  Var[R(k)  R(k-l)  ]  +  Var  [I  (k)  I  (k-1)  ] 

(3.18) 
Var{Im[Wpfk)  ]  }=  Var(Ip(k)}  =  Var  [I  (k)R(k-l)  ]  +  Var[R(k)  I  (k-1)  ]  . 

Let  X  and  Y  be  statistically  independent,  non-zero  mean  random 
variables.  Then 


17 


Var[XY]  =  E[xV]  -  {E[XY]}2 


E[X2]E[Y2]  -  {E[X]}2{E[Y]}2 


o*oy2   +{E[X]}2ay2  +  {E[Y]}2ax2.  (3.19) 


It  can  be  seen  that  the  variance  of  a  product  contains  a  noise 
times  noise  term  and  two  signal  times  noise  terms;  the  later  two 
are  due  to  the  mean  values  of  X  and  Y.  Recalling  now  that  the  real 
and  imaginary  values  R(k)  and  I(k)  have  been  normalized  to  unit 
variance  for  all  k,  then  it  follows  that 


a2  =  Var{Rp(k)}  =  2(1  +  Ek/N0  +  EM/H0) 


a2  =  Var{I_(k)}  =  2(1  +  Ek/N0  +  EM/N0) 


(3.20) 


Here,  the  first  term  is  due  to  noise  times  noise,  the  second  two 
are  the  signal  times  noise  terms.  This  is  important  because  at  high 
SNR,  the  signal  times  noise  terms  dominate  and  R  amd  I  are 
approximately  gaussian  with  mean  values  given  by  (3.16)  and 
identical  variances  given  by  (3.20). 

It  is  important  to  determine  if  the  real  and  imaginary  parts 
of  the  complex  product  noise  terms,  W  (k)  are  correlated.  We  note 
that,  because  W  (k)  is  zero  mean, 

Covar{RD(k)l(k)}   =  E{Re[WD(k)  ]  Im[WD(k)  ]  } .     (3.21) 
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The  covariance  is  defined  as 

Covar{Rp(k)Ip(k)}  =  E{Rp(k)  Ip(k)  }  -  E{iyk)  }E{Ip(k)  }  .  (3.22) 

which  we  find  by  using  (3.17)  and  the  previously  cited  properties 
of  the  noise  to  be  identically  zero (this  is  so  at  all  SNRs  and  is 
not  dependent  on  the  gaussian  approximation) .  Zero  covariance 
insures  that  D(k)  ,  the  rotated  complex  product  of  (3.13)  also  has 
uncorrelated  real  and  imaginary  parts  with  equal  variances  as  given 
by  (3.20).  At  high  SNRs  it  tends  to  the  gaussian. 

Dividing  (3.13)  by  a  produces  a  normalized  differentially 
demodulated  complex  random  variable  D' (k)  with  unit  variance  real 
and  imaginary  parts  and  mean  value  given  by, 


E[D'(k)]=  (S'(k)S'  (k-l)/ap}exp{jA0/2}  = 


M'kcos(d0(k)+A0/2)  +jM'ksin(d0(k) +A0/2)   (3.23) 
where 

M\  =[(2EkEk.1)VN0]/[l+Ek/N0+Ek.1/N0],/2.  (3.24) 

Comparing  (3.23)  and  (3.24),  with  (1.4)  and  (1.5)  it  can  be  seen 
that  differential  modulation  causes  a  reduction  in  SNR  for 
decoding.  If,  for  example,  adjacent  tones  have  equal  energy  and  if 
the  SNRyB  =  Ek/N0  is  much  larger  than  one,  then 
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M'k  *  <Ek/N0}*  =pk  (3.25) 


which  is  3  db  less  than  (1.5)  predicts  for  the  coherent  PSK. 

The  symbol  error  probabilites  are  determined  as  before.  So 
long  as  the  SNR  is  sufficiently  high  to  use  the  gaussian 
approximation,  then  the  probability  of  a  symbol  error  for  MFDQPSK 
is 


Ps  =  1  -  [1  -  Q{M'k/(2)*}]2   *  2Q{M'k/(2)V  (3.26) 


The  probability  that  D(k)  falls  outside  a  sector  of  width  A6 
is 

p(A9)  »  2Q(M'ksin(A9/2)  }  (3.27) 


M 

so  the  probability  of  a  differential  symbol  decoding  error  in  2  ■ 
ary  MFPSK  is 


Ps  =  P(A0)  (3.28) 


with  M'k  given  by  (3.24)  and  A0  =  27r/2M. 


Bit  Error  Rates  for  MFPSK  and  MFDPSK 

In  order  to  compute  BERs,  the  expected  number  of  bit  errors 
per  symbol  must  be  computed  and  divided  by  the  number  of  bits  per 
symbol  (see  (2.10)).  The  probability  of  single,  double,  triple  bit 
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errors  etc.  depends  on  the  transformation  that  is  used  ,  that  is  it 
depends  on  (3.1), (3. 10), 


ELJUU-i 


i(k) A0. 


(3.29) 


Because  the  most  probable  errors  are  in  the  sectors  that  are 
adjacent  to  the  sector  containing  the  transmitted  symbol,  and  the 
next  most  probable  are  in  the  next  set  of  adjacent  sectors,  etc., 
we  have  elected  to  use  a  reflex  or  Grey  code  to  code  the  phase 
bits.  In  this  way,  adjacent  sectors  only  differ  by  one  bit,  so 
that  the  most  probable  symbol  decoding  errors  are  single  bit 
errors,  next  most  adjacent  sectors  differ  by  two  bits,  so  they 
yield  double  bit  errors.  In  this  way,  we  hope  to  minimize  BER  when 
the  errors  are  due  to  AWGN. 

The  Grey  reflex  code  for  2,  3,  and  4  phase  bits  is  given  for 
reference  in  Table  3,1 


ELOU-i 


i(k) 


M=2 

M=3 

M=4 

M=2 

M= 

=3 

M=4 

00 

00,0 

00,00 
00,01 

0 

0 

0 

1 

00,1 

00,11 
00,10 

1 

2 
3 

01 

01,1 

01,10 
01,11 

1 

2 

4 
5 
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01,0 


11    11,0 


11,1 


10    10,1 


10,0 


01,01 
01,00 
11,00 
11,01 

11,11 

11,10 
10,10 
10,11 
10,01 
10,00 


6 

7 

8 

9 

10 

11 

12 

13 

14 

15 


TABLE  3.1  Reflex  Coded  Phase  Bits 

Note  from  Figure  2  that  the  two  most  significant  bits  are 
determined  soley  from  the  quadrant  in  which  the  received  complex 
point  occurs.  The  remaining  bits  depend  only  on  the  sector  occupied 
by  the  absolute  values  of  the  real  and  imaginary  parts  of  the 
complex  values.  These  observations  are  useful  for  designing 
encoding  and  decoding  algorithims  for  variable  numbers  of  bits  of 
phase  modulation  in  a  systematic  fashion  (see,  Mercury  Digital 
Communications  patent  application  [7]). 

a.  MFQPSK  and  MFDQPSK  Bit  Error  Rates. 

It  is  easy  to  see  that  the  two  quad  bits  are  decoded  as 
follows; 
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Most  Significant  Bit  -»  o  if  I(k)  >  0; 
Most  Significant  Bit  ■*  1  if  I(k)  <  0. 

Next  Most  Significant  Bit  -»  0  if  R(k)  >  0; 
Next  Most  Significant  Bit  ■»  1  if  R(k)  <   0. 

Here  R(k)  and  I(k)  are  the  real  and  imaginary  parts  of  Y(k)  or  of 
D(k)  for  coherent  MFQPSK  or  MFDQPSK  respectively.  Because,  in  both 
cases  they  are  statistically  independent  random  variables,  then 
QPSK  is  equivalent  to  sending  two  independent  one  bit  symbols  and 
consequently  the  BER  is  just  the  symbol (bit)  error  probability  of 
either  one.  Thus,  for  QPSK 

BER  =  Q(X)  (3.30) 


where 


for  Coherent  MFQPSK; 


x  =  K\/(2)h     =  p  (3.31) 


and  for  MFDQPSK; 


x  =  M'k/(2)*  =  P/2*  (3.32) 
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b.  2M-ary  Coherent  MFPSK  and  MFDPSK. 


Due  to  the  symmetry  of  the  phase  sectors  and  the  Grey 
coding  algorithim  employed  in  the  PSK  and  DPSK,  it  is  clear  that 
the  conditional  probabilities  of  B  bit  errors  in  a  symbol  of  M 
phase  bits  are  all  equal.  That  is,  Pr[B=b|s  =s,-]  is  the  same  for 
all  s,- .  and  therefore  so  is  E[B|s  =  sj.  Thus,  for  PSK  systems 

BER  =  E[B|s]/M  (3.33) 


and 


E[B|s]  =  l{Pr[B=l|s])  +  2{Pr[B=2|s]}  +  ...+  M{Pr [B=M| s] } .  (3.34) 

Figure  3  illustrates  B  ,  the  number  of  bit  errors  for  M  =4  in 
each  of  the  16  sectors  with  respect  to  the  transmitted  symbol  s  = 
00,00.  Note  that  of  the  4  single  bit  errors,  2  are  in  adjacent 
sectors,  and  of  the  6  double  bit  errors,  2  are  in  sectors  adjacent 
to  the  single  bit  sectors.  We  shall  assume  for  the  remainder  of 
this  report  that  the  probability  a  received  complex  value  is  more 
than  two  sectors  away  from  the  transmit  sector  is  negligible 
compared  with  the  probability  of  being  within  the  two  closest 
sectors.  Thus, 

Pr[B=l|s]  «  p(A0)  -  p(3A0) 

(3.35) 
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Pr[B=2|s]  «  p(3A0) 


and  consequently 


BER  «{p(A0)  +  p(3A0)}/M, 


(3.36) 


with  p(A6)  determined  by  (3.8)  for  coherent  MFPSK  and  by  (3.27)  for 
MFDQPSK. 


Equations  (3.30) 
and  (3.36)  represent  the 
principal  results  of 
Chapter  HI.  They  are 
used  to  compute  BER  for 
coherent  phase  coded  and 
for  differentially  phase 
coded  MFM.  Grey  coding 
of  the  phase  bits  into 


i 

4 

\    \  1 

A"\ 

the  phase  sectors  has  Fi(?ure  3  Distribution  of  Bit  Errors  B  for 
*  s  =00,00. 


been  assumed. 
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IV  AMPLITUDE  MODULATION 

In  this  Chapter,  results  are  presented  for  digitally  coding 
information  into  the  amplitude  of  MFM  carrier  tones.  For  the  non- 
differential  case,  these  are  but  a  recitation  of  well  known  results 
for  incoherent  amplitude  shift  keying  since  all  the  tones  are 
orthogonal  and  are  in  effect  processed  in  parallel.  The  advantage 
of  amplitude  modulation  is  that  amplitude  coding  does  not  require 
any  knowledge  of  the  individual  carrier  tone  phases  for 
demodulation.  It  does  require  knowledge  or  control  of  the  gain  of 
the  channel  at  each  of  the  carrier  tones  in  order  to  set  the 
thresholds  for  the  decision  regions.  Consequently,  a  technique  for 
differentially  encoding  data  as  the  change  in  amplitude  between 
adjacent  carrier  tones  will  also  be  described  and  analyzed.  This 
method  and  an  apparatus  for  its  implementation  have  been  described 
in  a  US  Patent  Application[7]  .  So  long  as  the  channel  gain  is 
nearly  identical  at  adjacent  frequencies,  the  differential 
amplitude  coding  system  is  independent  of  both  the  phase  and  the 
amplitude  response  of  the  channel. 

Multilevel  Amplitude  Modulation 

Multilevel,  multi-frequency  amplitude  modulation  encodes  N 
amplitude  bits  into  one  of  2  amplitude  levels  in  each  of  the 
carrier  tones.  The  coding  is  prescribed  by  a  transformation- 


al  -  Ak  =  i(k)AA  +  A^  ;  i  =  0,1...2N-1,   (4.1) 
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and  as  with  the  MFPSK  coding,  a  Grey  code  is  used  so  that  adjacent 
amplitude  values  correspond  to  only  one  bit  difference  in 
codewords . 

Given  the  transmission  of  codeword  a(k) .  the  normalized 
decoded  complex  DFT  coefficient  F(k)  is  a  unit  variance  gaussian 
complex  random  variable  with  uncorrelated  real  and  imaginary  parts 
and  average  values  given  by  (1.5).  The  magnitude  of  F(k) ,  is  a 
Rice-Nakagami  random  variable,  R,  with  probability  density  function 


fR(r|a(k))  =  rexp{-(r2+A'k2)/2}I0(rA«k)         (4.2) 


Now,  for  A'k  >>1,  integrals  of  fR(r)  in  the  vicinity  of  A'k  can  be 
approximated  by  integrals  under  the  gaussian  distribution  with  mean 
value  A'k  and  unit  variance,  (see  e.g.  [9],  pg  504  and  pg  625.) 

Amplitude  bits  will  be  decoded  in  accordance  with  the  MLD  rule 
as  illustrated  in  Figure  4.  Conditional  amplitude  symbol  error 
probabilities  are; 


PS[A\  =  A'Bin]  *  Q(AA'/2) 


Ps[A'min    <A'k    <   A'MX]    «    2Q(AA'/2)  (4.3) 


Ps[A\   =   A'™  3    *   Q(AA'/2) 


Assuming  equally   likely  symbols,   the  average   symbol   errror 
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probability   is 


Ps   =    2"N   S    P8[A'k]    =2(1-    2"N)Q(AA'/2)  (4.4) 
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Figure  4  Multi-Level  Amplitude  Decoding 


Average  SNR^.E/Np  and  Eb/N0 

When  amplitude  modulation  is  used,  the  carrier  tones  have 
different  power.  However,  we  would  like  to  evaluate  the  error 
probabilites  in  terms  of  an  average  SNP^g  and  average  Eb/N0. 
Let  A'0  =E[Ak']  be  the  average  signal  amplitude  of  tone  k  and  let 
the  increment  between  amplitudes  be  proportional  to  the  average 
amplitude,  that  is  let  AA'=aA'0.  Then  it  is  easily  shown  that  the 
average  tone  energy  E  is  given  by 


2  ,„N 


E  =E[Ek]  =  E0{1  +(a/2)fc(2  -1)  [l+(2-2)/3]}    (4.5) 
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where 


E0  =  (A'0)2N0/2  (4.6) 


is  the  energy  associated  with  the  average  amplitude  of  tone  k.  We 
shall  be  most  interested  in  zero,  one,  two  and  three  bits  of 
amplitude  coding  for  which  N  is  zero,  one,  two  and  three 
respectively  and  (4.5)  shows  that  the  average  energy  is  related  to 
the  energy  of  the  average  amplitude  tone  by 


N=0:  E  =  E0 


N=l:  E  =  EQ{1  +(a/2)2} 


(4.7) 


N=2:  E  =  E0{1  +5(a/2)2} 


N=3;  E  =  E0{1  +21(a/2)2) 


The  average  SNRyB  and  average  E^Ng  are  related  to  the  average  tone 
energy  E  as  in  (1.6)  and  (1.7),  viz; 

SNRyg  =  E/No  "  P  (4.8) 


and 


VNo  =  E/QNo-  (4-9) 
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In  order  to  evaluate  symbol  error  probabilites,  we  need  to 
evaluate  Q(AA'/2).  Now  recalling  that  AA'  has  been  defined  as  a 
fraction  of  the  average  tone  amplitude,  then  it  follows  that 


Q(AA'/2)=Q(aAV2)  =  Q(a{2E0/N0)V2) 


=  Q({2E/N0}*f(a,N))  (4.10) 


with 


f(a,N)  =  a/[2{l+(a/2)2(2N-l)[l+(2N-2)/3]}*]     (4.11) 


It  is  clear  that  the  function  f (a,N)  is  of  central  importance 
for  controlling  error  probabilites.  Notice  that  there  is  a  maximum 
value  for  a  because  the  smallest  amplitude  level  must  be  non- 
negative  and  that  the  maximum  value  decreases  with  increasing  N. 
That  is, 


0  <  a  <  2/(2N-l)  .  (4.12) 


We  note  that  f(a,N)  can  be  reasonably  approximated  by  a/2  for 
values  of  a  less  than  one.  It  is  also  clear  that  if  one  is  just 
using  amplitude  modulation,  a  should  be  made  as  large  as  possible, 
which  is  the  same  as  spreading  the  tone  amplitude  levels  as  much  as 
possible,  in  order  to  minimize  error  probabilities.  Consequently, 
with  amplitude  modulation  alone,  performance  is  limited  by 


f(amax'N)  =  {l/(2N-l)}{l/[l+(l/(2N-l))(l+(2N-2)/3)]*}   (4.13) 
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For  one,  two,  three  and  four  bits  of  amplitude  modulation,  (4.13) 
is  given  in  Table  4.2  below. 

N     flQW^Nl 


1 

1/2* 

2 

1/14* 

3 

1/70* 

4 

1/310* 

TABLE  4.2 
The  Maximum  of  the  Spread  Function 

Finally,   combining   (4.4)   and   (4.10),   the   symbol   error 
probability  for  amplitude  modualtion  alone  of  the  tones  is  given  by 


P,  =  2(l-2"N)Q{(2)*pf(a,N)}.  (4.14) 


Amplitude  and  phase  modulation  symbol  error  probabilites  for 
different  numbers  of  bits  may  be  compared  by  comparing  (4.14)  with 
(3.9).  (  It  should  be  noted  that  for  the  situation  of  a  =  cnmax,  that 
the  minimum  amplitude  A'min=0  and  therefore  the  probability 
distribution  is  Rayliegh,  not  Gaussian  as  assumed  in  (4.3).  Since 
we  are  primarily  interested  in  combined  amplitude  and  phase 
modulation  constellations  for  MFM  for  which  0  amplitude  tones  are 
not  permissible,  we  shall  not  be  concerned  with  presenting  the 
necessary  correction  here.) 
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BER's  for  Amplitude  Modulation 

Bit  error  rates  for  the  amplitude  bits  are  determined  as 
presented  in  Chapter  II.  That  is, 

BER  =  E[B]/N  (4.15) 

with 

E[B]  =  Z  b  {Z  Pr[B=b|a]Pr[a] }   (4.16) 

We  shall  approximate  (4.16)  using  single  and  double  bit 
errors.  This  is  a  good  approximation  by  virtue  of  the  Grey  coding 
of  the  amplitudes  bits.  Thus,  assuming  equally  likely  amplitude 
codewords, 


E[B]  =  2"N  {S  Pr[B=l|a]  +  2  S  Pr[B=2|a]}    (4.16) 


These  conditional  error  probabilities  depend  on  the  whether,  the 
codeword  a  corresponds  to  A'min  or  A'max  or  to  an  in-between  value 
of  A1 .  We  are  principally  concerned  with  one,  two  and  three  bits 
of  amplitude  modulation.  The  evaluation  of  (4.16)  for  these  cases 
is  the  following: 

N=l: 

E[B]  =  Ps  =  Q(AA'/2)  =  BER  (4.17) 

N=2: 

E[B]  =  1.5Q(AA'/2)  +  Q(3AA'/2)=  2BER  (4.18) 
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N=3: 

E[B]  =  (7/4)Q(AA'/2)  +  (3/2 ) Q (3AA'/2)  =  3BER   (4.19) 


with 


AA'/2  =  (2}*pf(a,N)  (4.20) 


Differential  Amplitude  Modulation 

In  order  to  decode  the  amplitude  bits  of  MFM  carrier  tones,  it 
is  necessary  to  establish  accurate  values  for  the  channel  gain  at 
each  carrier  tone  frequency.  Although  this  is  less  difficult  to  do 
than  to  determine  phase  shift,  it  may  still  be  of  advantage  to  have 
a  modulation  scheme  which  does  not  require  either  phase  or 
amplitude  information  in  order  to  decode  the  symbols.  In  mobile 
systems  with  bandpass  signals  subject  to  fading  this  will  almost 
certainly  be  the  case.  Also,  in  switched  systems  where  a  premium  is 
on  connect  time,  allocation  of  valuable  connection  time  to 
measurement  of  channel  gain  may  represent  a  limitation. 
Consequently,  we  have  investigated  the  performance  of  MFM 
modulation  with  differentially  encoded  amplitude  information. 
Although  encoding  more  than  one  bit  differentially  is  conceptually 
possible,  to  date  we  have  only  analyzed  and  tested  modulation 
schemes  involving  one  bit  of  differential  amplitude  coding  where 
the  coding  is,  as  with  the  differential  phase  modulation, 
differential  between  adjacent  tones.  Thus,  two  amplitudes,  A^  and 
A^  that  are  separated  by  AA  =  aA0  =  a{\ax    +    Amin}/2  are  used. 
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Differential  encoding  the  amplitude  bit  of  tone  k  is  accomplished 
by  making  Ak  the  same  as  Ak.1  if  a=0  and  changing  Ak  to  the  amplitude 
different  from  the  amplitude  of  Ak.1  if  a=l.  The  decoding  rule  for 
a  is  illustrated  by  Figure  5.  The  crosshatched  region  is  decoded  as 
a=0,  that  is  it  is  assumed  that  the  amplitudes  are  the  same  in  this 
region.  The  remaining  region  is  decoded  as  a=l,  that  is  it  is 
assumed  that  the  amplitudes  differ  in  this  region.  This  rule  is 
described  by  the  algorithim; 

"if    (l/(l-a/2) } |Y(k-l) |<|Y(k) |<(l-a/2) |Y(k-l) |  then  a=0 
else  a=l" . 


Clearly  the  rule  is 
independent  of  channel  gain  so 
long  as  adjacent  tones  have 
essentially  identical  gain. 

To  find  symbol  and  bit 
error  probabilities,  the 
conditional  probability  of  |y| 
being  outside  the  decision 
regions  must  be  computed.  Recall 
that  we  have  established  the 


i 

>K 

/   / 

/  /       / 

0*JL_      /  / 

/ .         /          / 

■km      •Ay    /y 

' 

'/     /        jS 

//    /     js 

-Km  //    •    ^X      • 

///yr      ^ 

/     Ayy 

J^/ 

r    s          *.               * 

V 

kr*\*      **       ki**t 

kk-{ 

mean  and  variance  of  |  Y  |  and  the  Figure  5  Differential  Amplitude 

Decoding 

fact  that  adjacent  magnitudes 

are  statistically  independent.    Consider  the  dotted  lines  with 

slope  of  one  passing  thru  the  points  (A0,Am1n)  and  (ABin,A0)  .  We  shall 
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use  these  as  approximations  for  actual  decision  boundaries  shown  by 
the  solid  lines.  Accordingly,  it  can  be  seen  that  the  average 
symbol  error  probability  which  is  also  the  BER  is  given  by; 


BER  =  Ps  =  (3/2)Q{  (l/2*)AA'/2}  -  (1/2)  Q{  (l/2)*3AA'/2  }  .   (4.21) 


Comparing  (4.21)  and  (4.17),  we  note  that  a  differentially  encoded 
amplitude  bit  requires  approximately  3db  more  signal-to-noise  ratio 
than  an  absolutely  encoded  amplitude  bit  in  order  to  obtain  the 
same  BER. 
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V.  QUADRATURE  AMPLITUDE  MODULATION 

The  purpose  of  this  Chapter  is  to  describe  more  complex 
modulation  constellations,  constellations  that  combine  the 
amplitude  and  phase  modulation  schemes  discussed  in  the  previous 
two  Chapters.  As  mentioned  in  the  Introduction,  we  shall  refer  to 
these  as  "MFQAM  systems"  even  though  they  will  not  have  Z2  (that  is 
rectangular)  grid  points  as  is  commonly  found  in  single  carrier  QAM 
constellations . 

To  be  more  specific,  we  will  first  consider  systems  that 
utilize  differential  phase  modulation  of  M  bits  and  one  bit  of 
differential  amplitude  modulation.  These  systems  will  be  referred 
to  as  "fully  differential  systems".  Also,  we  will  consider  systems 
that  use  differential  phase  coding  and  absolute  amplitude  coding. 
A  system  of  this  type  with  M  bits  of  differential  phase  modulation 
and  N  bits  of  absolute  amplitude  modulation  will  be  referred  to  as 
anrt(M,N)  Hybrid". 

Single  Amplitude  Bit  MFQAM  Systems 

These  systems  accept  Q=M+1  bit  codewords  s  which  are 
partitioned  into  two  quad  bits  that  determine  the  phase  difference 
of  adjacent  tones  to  the  nearest  ninety  degrees  ,  M-2  delphase  bits 
that  are  coded  into  phase  within  the  quadrant  and  one  amplitude  bit 
the  is  encoded  as  tone  amplitude.  Thus  a  codeword  is  divided  into 
the  three  subsets 
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s  =  g,d,a  .  (5.1) 

As  described  in  Chapters  II  and  III,  all  coding  is  done  in 
accordance  with  a  Grey  (reflex)  code  such  that  nearest  neighbor 
points  in  QAM  constellations  differ  in  only  one  bit  position.  The 
BER  for  these  constellations  is  determined  from  (2.10)  using  (2.11) 
and  (2.12).  We  shall  consider  only  single  and  double  bit  errors. 
Accordingly, 

Pr[B=l|s=si]  and  Pr  [  B=2  |  s=s,.  ]  (5.2) 


are  to  be  evaluated  for  all  s,-. 


a. Single  Bit  Errors 

Single  bit  errors  occur  when  an  amplitude  bit  decoding  error 
occurs  and  all  the  phase  bits  are  decoded  correctly  or  when  a 
single  bit  phase  error  occurs  and  the  amplitude  bit  is  decoded 
correctly. 
These  correspond  to  the  error  patterns; 

e  =  0,1  and  e  =  ep1 ,  0  (5.2) 

where  e^1  is  a  single  bit  phase  error  pattern. 
Therefore  the  single  bit  error  probability  is, 

Pr[B=l]  =  Pr[ea=l,ep=0]  +  Pr[ep1,ea=0]     (5.3) 
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If  we  assume  amplitude  and  phase  decoding  errors  are  independent 
(which   is   certainly   not   strictly   true   but   is   probably   a 
satisfactory  approximation  at  low  error  rates)  then, 

Pr[B=l]  =  Pr[ea=l]Pr[ep=0]  +  Pr^]  { 1-Pr  [ea=l]  } 

*  Pr[ea=l]  +  Prfe^]  (5.4) 

where  the  approximation  is  accurate  at  low  symbol  error  rates. 

The  probability  of  an  amplitude  bit  error  is  clearly 
independent  of  the  phase  of  tones  and  is  given  by  (4.17)  and  (4.21) 
for  absolute  and  differential  amplitude  coding  respectively. 

The  probability  of  phase  bit  errors  is  identical  for  all 
phases  but  depends  on  the  amplitude  of  the  complex  rotated  product 
which  we  have  called  Mk  and  whose  statistics  are  given  in  (3.24). 
A  single  amplitude  bit,  that  is  two  amplitude  levels,  results  in 
three  different  complex  rotated  product  amplitude  levels  for  Mk 
depending  on  whether  adjacent  tones  both  have  amplitudes  A^,  both 
have  amplitudes  Amax,  or  one  has  amplitude  A^  and  the  other 
amplitude  Ainax.  Since  we  have  assumed  all  symbol  patterns  are 
egually  likely,  then  the  first  two  situations  each  occur  with 
probability  1/4  and  the  latter  with  probability  1/2. 

Let  us  designate  the  three  levels  of  M'k  as  M'^-^M'^  and  M'max. 
Then  the  probability  of  any  particular  phase  bit  error  pattern  is 

Pr[ep]  =  ^Pr[ep|M'mi.d]  +  ^<Pr  [ep|  M'min  +  Pr  [ep|M'max]  )   (5.5) 
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The  probability  of  single  (and  double)  phase  bit  error  patterns  are 
given  in  (3.35)  where  p(A0)  is  given  in  terms  of  M'k  by  (3.27)  for 
MFDPSK.  The  three  levels  for  M'k  are  found  from  (3.24)  for  E/N0  > 
1  to  be, 


"'in  =  (WN0>* 


"'id    =     (EminE.ax/ENo}"  (5.6) 


M'™    =    {Emax/N0}%. 


In  terms  of  the  amplitude  level  spread  parameter  a, 
Emin'A  =  Eo'^1"01/2)'  Emax%  =  E0*(l+a/2)  and  E0,  the  energy  of  an  average 
amplitude  tone  is  related  to  the  average  energy  of  the  tones  by 
(4.5).  Specifically  we  find  that, 


M'min  =  <(E/N0)  (l-a/2)/(l+a/2)>* 


M'mid  =  <(E/N0)  (l+(a/2)2)}*  (5.7) 


M'^  =  {(E/N0)(l+a/2)/(l-o/2)}*. 


Thus  (5.7)  in  conjuction  with  (5.5)  and  (3.35)  determines  the 
probability  of  single  phase  bit  errors  in  terms  of  the  p=  SNRyB  = 
E/N0,  the  number  of  phases  2M  and  the  parameter  a.    Finally  the  BER 
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due  to  single  bit  errors  is  well  approximated  by 


BER,  =Q"1{Pr[ea=l]   +  Prtfip,])  (5-8) 


with 


Q(pa0)     for  hybrid  coding 


Pr[e=l]  = 


1(3/2)  Q(pa0/2%)    -(l/2)Q(3pa0/2*)    for   fully 


differential  coding,  and 


Prtep,]   =  Q(pa2)-Q(3pa2)+^[Q(pa1)-Q(3pa1)+Q(pa3)-Q(3pa3)] 


with 


p   =  (E/NJ* 


aQ  =  a/<2[l+(a/2)2]}* 


a,    =   sin(7r/2M){(l-a/2)/(l+a/2)}V? 


a2   =   sin(7r/2M)  {l+(a/2)2}% 
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a3  =  sin(7r/2M)  {  (l+a/2)/ (l-a/2)  }*" 


Examination  of  (5.8)  shows  that  all  the  a  values  increase  with  a 
except  a1  which  determines  the  phase  errors  due  to  complex  products 
of  the  minimum  amplitude  tones.  This  makes  sense  as  we  know  that 
the  amplitude  errors  decrease  with  increasing  amplitude  difference 
and  the  large  amplitude  tones  will  have  more  energy  with  increasing 
amplitude  differential. 

One  design  procedure  we  have  investigated  is  to  choose  a  such 
that  a1  =  a0  for  the  hybid,  i.e.  absolute  amplitude  bit  coding  and 
a1  =  a0/2  for  differential  amplitude  coding.  This  leads  to  the 
following  values  of  a; 


M 


a   Fully  Differential  Coding 


a   Hybid  Coding  a   Geo 


3  2 

4  3 

5  4 

6  5 


.938 
.593 
.335 
.180 


728 
447 
250 
130 


.88 
.54 
.32 
.18 


TABLE  5.1 
Amplitude  Spread  Factor  for 
a  Single  Bit  Encoded  as  Amplitude 


From   a   purely   geometrical   point   of   view,   values   of 
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a=A0/ (1+A0/2)  space  the  constellation  points  approximately 
equally.  These  values  are  listed  for  comparison  in  Table  5.1  too. 
Using  the  optimum  values  of  a  yeilds  the  a  coefficients  shown  in 
Table  5.2. 


M=2 

Differential 

.601 

.425 

.778 

1.20 

Hybrid 

.485 

.485 

.740 

1.05 

M=3 

Differential     .403       .285      .410 

Hybrid  .310       .310       .40 


.525 
.482 


M=4 

Differential    .232 

Hybrid  .172 


.164 
.172 


.20 
.20 


.231 
.217 


M=5 

Differential    .127 

Hybrid         .092 


.090 
.092 


.10 
.10 


.110 
.107 


TABLE  5.2 
Error  Function  Coefficients 
for  MF-QAM 
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Theoretical  values  of  SNRyB   =  E/N0  =  p   are  given  in  db  in 

-2       -4-6  . 

Table  5.3  for  BERs  of  10   ,  10   and  10   due  to  single  bit  errors 
for  MF16-QAM.  (E^Nq  values  are  6  db  lower.) 


10"2       10"*       10"6 


Differential        17.9      21.9      24.2 

Hybrid  15.9      21.0      23.3 

TABLE  5.3 
Theoretical  SNRs  in  db  for  BER^  shown 
for  Fully  Differential  and  Hybrid  MF16-QAM 


b.  Double  Bit  Errors 

Double  bit  errors  are  less  likely  than  single  bit  errors  but 
can  contribute  to  the  overall  BER  at  the  lower  values  of  SNRyB. 
The  probability  of  a  double  bit  error  is 

P[B=2]  =  Pr[ep2/ea=0]  +  Pr[ep1,ea=l] 

<  Pr[ep2]  +  Pr[ea=l]  (5.9) 
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where  the  last  term  in  the  inequality  is  used  as  a  bound  because 
of  the  probable  high  correlation  between  amplitude  and  phase  bit 
errors  and  is  given  by  (5.8) .  The  double  phase  error  probability 
is  determined  from  (5.5)  and  (3.35)  with  (3.27).  This  is  always 
much  less  than  the  probability  of  an  amplitude  bit  error  when  the 
amplitude  spread  has  been  adjusted  to  a  value  such  that  amplitude 
bit  errors  and  single  phase  bit  errors  are  approximately  equal. 
Thus  we  may  approximate  the  BER  due  to  double  bit  errors  as 


BER2  «  2Q"1Pr[ea=l]  (5.10) 


with 


Q(pa0)     for  hybrid  coding 


Pr[e=l]  = 


;(3/2)Q(pa0/2Vj)    -(l/2)Q(3pa0/2%)    for    fully 
differential  coding. 

The  effect  of  double  bit  errors  is  to  approximately  double  the 
BER  determined  by  BER1 .  However,  because  the  Q(x)  function  falls 
off  so  rapidly,  this  can  be  offset  by  0.3,  0.5,  and  0.8  db 
increases  in  SNR  for  BERs  of  10  ,10"  ,  and  10   respectively. 

Multi-level  Amplitude  MFOAM  Systems 

These  systems  become  quite  complex  to  analyze  exactly  because 
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of  the  multiplicity  of  amplitudes  that  occur  in  the  complex  rotated 
product  demodulation  of  the  differential  phase  terms.  However,  a 
reasonable  hueristic  design  procedure  is  to  set  the  differential 
amplitude  level  such  that  the  distance  in  the  constellation  between 
the  two  smallest  amplitudes  is  equal  to  the  radial  distance  between 
adjacent  phase  points  on  the  smallest  amplitude  ring.  This  is 
illustrated  in  Figure  6.  The  design  procedure  is  to  set 

AA'  =  A'min  A0  (5.11) 

(which  leads  to  a  ={1/A0  +  (2N-l)/2  }'1)  . 

With  this  constraint,  we  see  that  the  single  bit  error  probability 

will  be  bounded  by 

P[B=1]  <  Pr[ea=l]  +  Pr[epl]         (5.12) 


which  can  be  found  from  (4.3)  for  the  amplitude  bit  errors  and  from 
(3.35)  and  (3.27)  for  the  phase  bit  errors.  We  now  use  in  addition 
the  facts   that 


AA'  =  2  (2E/N0}%  f(a,N)  (5.13) 


with   f(a,N)    given  by    (4.11)    and   A0  =   27r/2M.    Thus, 


Pr[e=l]    =    2Q(AA'/2) 
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and 


Prtfip,    ]    =   2Q(M'm,n   A0/2) 


with  M'min  =    (Emin/N0)%  =  A'min/2%  =   AA'/2*A0   such   that 


PrCgp,]    =    2Q{AA'/(2(2  ))} 


(5.14) 


It  can  be  noted  that  the  phase  errors  are  3  db  worse  than  the 
amplitude  errors  because  of  the  differential  coding.  However,  the 
amplitude  errors  are  the  same  for  all  amplitude  levels  whereas  the 
phase  errors  are  much  less  for  the  larger  amplitude  values.  Using 
this  procedure  one  obtains  the  parameters  of  Table  5.4  for  the 
(M,N)  hybrids  indicated. 


Figure  6 

Ring  Spacing  to  Equalize  Amplitude  and  Phase  Bit  Errors 
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(M.N)         a  f  (a.N)    =   AA'/(2(2*)o)     SNR  in  db  (BER) 

io"6.io"4.io  2 


28  24.8 

122  32.0 

075  36.3  ,34.3,31.0 

058  38.5 

046  40.5 


TABLE  5.4 
SNRs  Required  for  Single  BER 
shown  for  Hybrids  from  4  to  8  Bits/Hz. 


(3,1) 

.563 

(4,2) 

.248 

(5,2) 

.152 

(5,3) 

.116 

(6,2) 

.085 

The  (3,1)  16-QAM  case  has  been  included  for  comparison.  Notice 
how  these  geometrical  arguments  lead  to  results  very  close  to  those 
obtained  in  the  previous  section  by  more  thourough  analysis.  (24.8 
vice  23.3  db  for  10"6  BER  and  a  of  .563  vice  .447).  This  gives  us 
considerable  confidence  in  our  other  results.  The  (5,2)  hybrid 
BERs  also  seem  to  be  closely  verified  by  our  current  experiments. 
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VI  DISCUSSION  AND  CONCLUSIONS 

The  principal  objective  of  this  report  was  to  document  the 
theortical  predictions  for  various  types  of  differential  coding  in 
the  frequency  domain  for  Multi-Frequency  Modulation (MFM) .  In  order 
to  do  this  in  a  systematic  fashion,  general  results  were  derived 
for  probability  of  symbol  error  and  bit  error  rate(BER)  in  Chapter 
II  in  terms  of  e,  the  error  pattern  vector.  The  coding  we  have 
developed  partitions  a  symbol  s  into  phase  bits  p_  and  amplitude 
bits  a.  Phase  bits  are  further  partitioned  into  two  quad  bits  g 
that  determine  the  phase  quadrant  and  any  remaining  phase  bits  that 
determine  tone  phase  within  a  quadrant. 

In  Chapter  III  results  were  presented  for  coherent  phase 
modulation  with  maximum  likelihood  decoding  and  for  differential 
phase  modulation  using  the  "complex  product  rotation"  algorithim 
followed  by  maximum  likelihood  decoding.  In  order  to  minimize  BER, 
the  phase  states  are  Grey  coded  so  that  the  most  probable  symbol 
errors  result  only  in  single  bit  errors. 

Differential  phase  coding  and  decoding  using  complex  rotated 
product  results  in  a  theoretical  3  db  loss  of  performance  versus 
SNR  (or  E^Nq)  when  considering  detection  on  a  tone  by  tone  basis. 
A  well  known  phenomena  in  differential  coding  is  the  increased 
probability  of  error  pairs[8].  This  occurs  because  the  noise 
present  in  a  single  tone  affects  the  complex  rotated  product  for 
two  adjacent  symbols.  Thus  if  noise  is  high  in  one  tone,  there  is 
an  increased  probability  of  two  adjacent  symbol  decoding  errors.  We 
have  not  presented  a  theoretical  treatment  of  this  phenomena  in 
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this  report  but  will  in  a  subsequent  report  that  treats 
interleaving  and  forward  error  control  techniques,  including 
trellis  coding,  for  MFDM.  However,  it  is  mentioned  here  because  in 
addition  to  the  3  db  of  performance  loss,  error  pairing  is  another 
penalty  paid  for  differential  coding. 

In  Chapter  IV  results  were  presented  for  amplitude  coding  the 
tones  of  MFM.  Multi-level  amplitudes  with  equally  spaced  intervals 
were  assumed.  The  performance  of  maximum  likelihood  detection  of 
the  amplitude  interval  of  each  tone  at  the  receiver  was  given.  This 
technique  is  of  course  incoherent  because  it  does  not  require 
reference  phases  for  the  tones  for  decoding,  but  it  does  require 
knowledge  of  the  absolute  gain  of  the  channel  at  each  frequency  in 
order  to  set  the  decoding  thresholds.  Consequently,  we  also 
presented  results  for  differential  coding/detection  of  a  single 
amplitude  bit.  Again,  about  a  3  db  penalty  is  paid  for  differential 
coding  and  the  increased  probability  of  error  pairing  can  be 
expected  here  as  well. 

In  Chapter  V  the  concepts  of  phase  and  amplitude  coding  were 
combined  to  create  MFQAM  constellations.  All  of  these 
constellations  are  circularily  symmetric  and  because  their  phase  is 
differentially  encoded,  the  decoding  is  independent  of  absolute 
tone  phase  shift  thru  the  channel.  Constellations  with  a  single 
differentially  encoded  amplitude  bit  and  differential  phase  are 
termed  "fully  differential"  constellations.  Constellations  with 
absolute  level  decoding  of  amplitude  bits  but  differential  decoding 
of  the  phase  bits  are  termed  "hybrid"  constellations.  SNR  required 
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to  achieve  BER  of  10  were  given  for  a  number  of  the  constellations 
in  Table  5.4.  Constellation  design  procedures  were  presented  that 
purport  to  very  nearly  minimize  BER. 

The  SNRs  (or  E^N^  required  for  differential  constellations 
to  achieve  given  levels  of  BER  are  about  3  to  5  db  higher  than  is 
required  for  fully  coherent  MFM  using  rectangular  (Z2)  grid  for 
locating  the  tone  real  and  imaginary  parts.  The  advantage  gained  by 
differential  coding  is  the  removal  of  any  requirement  to  equalize 
the  channel  and  to  determine  reference  phases  for  each  of  the 
tones.  In  certain  applications,  the  penalities  of  differential 
coding  in  SNR  will  be  outweighed  by  its  simplicity  of 
implementation,  synchronization  and  robustness  in  the  face  of 
changing  channel  parameters.  In  other  applications,  it  will  not  be 
worth  the  penalty. 

The  purpose  of  this  report  was  to  present  design  procedures 
and  performance  expectations  for  frequency  domain  differentially 
encoded  MFM  with  signalling  efficiencies  of  2  to  8  bits/Hz. 
One  of  the  newest  techniques  of  error  control  in  single  carrier 
digital  modulation  is  trellis  coded  modulation (TCM)  [10]  .  Our  direct 
funded  research  project  is  now  directed  toward  applying  TCM  to  MFM 
by  trellis  coding  in  the  frequency  domain.  In  brief,  we  believe 
that  TCM  has  the  potential  for  improving  our  raw  Eb/U0  performance 
by  3  to  4  db.  Furthermore,  with  interleaving  it  is  hoped  that  the 
effects  of  differentially  induced  error  pairs  can  be  removed. 

Application  of  TCM  to  pure  phase  coding  is  believed  to  be 
relatively  straighforward.  From  an  implementaion  point  of  view,  the 
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MFM  receiver  must  be  equipped  with  a  Viterbi  decoder.  Although 
complex,  the  procedures  are  well  established  to  do  this  and  Viterbi 
decoders  for  TCM  have  been  written  for  the  Motorola  56000  DSP  chip 
[11].  Trellis  coding  the  QAM  constellations  analyzed  in  Chapter  V 
may  not  be  a  straightforward  application  of  known  results.  Previous 
work  has  concentrated  on  TCM  for  Z2  grids  [10]  and  the  conversion  to 
our  circularly  symmetric  constellations  is  not  directly  apparent. 
This  problem  is  at  the  top  of  our  current  research  agenda. 


51 


REFERENCES 

1.  B.  Hirosaki,  "A  19.2  kbps  Voice-band  Data  MODEM  Based  on 
Orthogonally  Multiplexed  QAM  Techniques",  IEEE  Int'l  Conf.  Comm. 
Rec,  pp  661-665,  Aug  1985. 

2.  B.  Hirosaki  et.  al.,  "Advanced  Group-band  MODEM  Using 
Orthogonally  Multiplexed  QAM  Techniques",  IEEE  Trans.  Comm.,  Vol 
Com-34,  pp  587-592,  June  1986. 

3.  L.  J.  Cimini,  Jr.,  "Analysis  and  Simulation  of  a  Digital  Mobile 
Channel  Using  Orthogonal  Frequency  Division  Multiplexing",  IEEE 
Trans,  on  Comm.,  Vol  Com-33,  No.  7,  July  1985. 

4.  M.  Alard  and  R.  Halbert,  "Principles  of  Modulation  and  Channel 
Coding  for  Digital  Broadcasting  for  Mobile  Receivers",  EBU  Review, 
No.  224,  August  1987. 

5.  T.  K.  Gantenbein,  Implementation  of  Multi-Frequency  Modulation 
on  an  Industry  Standard  Computer,  Naval  Postgraduate  School, 
Monterey,  CA. ,  MSEE(Sept  1989),  March  1990. 

6.  P.  H.  Moose,  "Theory  of  Multi-Frequency  Modulation  Digital 
Communications",  NPS  Tech.  Rpt.  NPS  62-89-019,  May  1989. 

7.  P.  H.  Moose,  "Multi-Frequency  Differentially  Encoded  Digital 
Communications  for  High  Data  Rate  Transmission  Through  Unequal ized 

52 


Channels",  U.  S.  Patent  Pending,  March  1990. 

8.  Simon  Haykin,  Digital  Communications,  John  Wiley  and  Sons,  1988. 

9.  Leon  Couch  11,  Digital  and  Analog  Communications  Systems,  3rd 
Edition,  MacMillan,  1990. 

10.  G.  Ungerboeck,  "Trellis-Coded  Modulation  with  Redundant  Signal 
Sets",  IEEE  Comm.  Magazine,  Feb.  1987. 

11.  Dion  Messer,  "Convolutional  Encoding  and  Viterbi  Decoding  Using 
the  DSP56001  with  a  V.32  Modem  Trellis  Example",  APR6/D,  Motorola 
Inc.,  Phoenix,  AZ . ,  1989. 


53 


DISTRIBUTION  LIST 


No.  Copies 


Defense  Technical  Information  Center  2 

Cameron  Station 
Alexandria,  VA  22304-6145 

Library,  Code  0142  2 

Naval  Postgraduate  School 
Monterey,  CA  93943-5002 

Department  Chairman,  Code  62  1 

Naval  Postgraduate  School 
Monterey,  CA  93943-5004 

Director  of  Research  Administration,  Code  012     1 
Naval  Postgraduate  School 
Monterey,  CA  93943-5000 

Paul  H.  Moose,  Assoc.  Prof,  Code  62me  12 

Naval  Postgraduate  School 
Monterey,  CA  93943-5004 

Darrell  Marsh,  Code  624  2 

Naval  Ocean  Systems  Center 
San  Diego,  CA  92152 


DUDLEY  KNOX  LIBRARY 


3  2768  00392858  1 


